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1> ■ Abstract 

We use singular Euclidean solutions to find multi-particle production cross 
^ ■ sections in field theories. We investigate a family of time-dependent 0(3) sym- 

' metrical solutions of the Yang-Mills equations, which govern the isotropic high- 

energy gauge boson production. At low energies our approach reproduces the 
instanton-induced cross sections. For higher energies we get new results. In 
particular, we show that the cross section for isotropic multiparticle produc- 
tion remains exponentially small in the running gauge coupling constant. The 
result applies both to the baryon number violation in the electro-weak theory 
and to the QCD jet production. We find that the isotropic multi-gluon pro- 
duction cross section falls off approximately as a ninth power of energy but 
possibly might be observable 



^Alexander von Humboldt Forschungspreistrager 



1 Introduction 



The last few years have witnessed an increasing interest in the non-perturbative 
multi-particle production induced by classical solutions of the field equations. 
Typical examples are (i) baryon number violation in the electro-weak theory 
(for reviews see ||l], (ii) multi-jet production in strong interactions ^ 

and (iii) multi-pion production in heavy ion collisions 0, 

While there seems to be a consensus in that classical field configurations 



are relevant in many many processes (see, e.g. ref. IC]), it is not 

so clear in the case of the 2 many ones. The difficulty with the latter 

processes is that an initial state with a few number of particles is a quantum 
rather than a classical one - even at asymptotically high energies. If one starts 
from a classical field with large energy at t = +cxd and evolves it back in 
time according to the classical equations of motion which preserve the field 
energy, one would end up with another field configuration at t = — oo with 
the same energy. Meanwhile, at t = — oo the whole energy has to belong 
not to the classical field (which should be zero or, at best, contain only the 
wrong-frequency part not corresponding to any real particles) but to the few- 
particle quantum state. Therefore there cannot exist any continuous classical 
field configuration interpolating in time between t = ±00. 



A way to overcome this difficulty has been indicated by Khlebnikov |11] 
who, following a remarkable work of lordansky and Pitaevsky |12], has sug- 
gested to study singular Euclidean solutions of field equations. Singular fields 
do not conserve energy across singularities, therefore they are adequate to de- 
scribe transitions between a quantum state with low or zero field energy and 
a (semi) classical one with high energy. Singular trajectories in imaginary (Eu- 
clidean) time have been introduced 60 years ago by Landau to calculate matrix 



elements between low- and high-energy states in quantum mechanics [13|. Even 
if a final high-energy state can be described semiclassicaly, the initial low-energy 
state needs not; therefore the Landau theory gives an example how, neverthe- 
less, the matrix elements can be treated semiclassically. It should be mentioned 
that an application of the Landau theory to multi-particle production has been 
suggested by Voloshin |^], however only constant fields and hence relatively 
low energies have been considered in that work. 

The energy of the initial state, i.e. the collision energy, is introduced as 
follows: Let the multiparticle production be initiated by, say, an annihilation 
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of two high-energy particles with energies £^1^2 in the cm. frame. To get the 
physical on-shell amplitude one has to calculate the two-point Green function in 
the singular background field, take its Fourier transform and then apply the leg 
amputation procedure of Lehmann-Symanczik-Zimmermann. The high-energy 
asymptotics of the Green function is apparently given by the Fourier trans- 



form of the singularity point [12|, which is exp {—i{Ei + E2)tsing) in Minkowski 
space, or, if one passes to the Euclidean space with the usual substitution 
—it — > t, one gets exp{Etsing) where E = Ei + E2 \s the total cm. energy of 
the process and tging is the time position of the field singularity in Euclidean 
space. We will see below that integrating over the positions of the singularity 
will result in the necessary conservation law: E = Efi^i^ where Efi^i^ is the 
energy of the produced multi-particle state. 

Our aim is to calculate semi-classically the total cross section induced by 
two high-energy particles with a total energy E, i.e. the imaginary part of a 
forward scattering amplitude which, in its turn, can be expressed through the 
4-point Green function (we use (f){r,t) to describe a generic field): 

a{E^ = (pi + p2f) = (flux factor) lim Im / d^xi_4e~'P'<'^^~''-^^~'P''-''''^~''^'> 

{df + m'^)(j)ixi){di + m^)(j){x2){dl + m^)(j)ix3){dl + m^)(l){x4)) . (1.1) 

The purpose of the paper is to demonstrate how to calculate the multi- 
particle production cross sections given generically by eq. (^) , using the tech- 
nique of singular Euclidean solutions. 



2 Singular Trajectories in Quantum Mechan- 



ics 



Before proceeding to the Yang-Mills case we describe a quantum mechanical 
analogue of our approach. We consider a particle with coordinate (j){t) in a 
double-well potential, whose (Euclidean) action is given by 



1 

7 



dt 



1 



1 



2^^ + 8< 
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(2.11 



where 5^ <C 1 is a dimensionless coupling constant. 
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2.1 Quantum-Mechanical "Cross Section" 



In the quantum- mechanical case " mass sheh" means fixed energy for each " par- 
ticle" , equal to 1 in our notations. Therefore, to get an analogue of high-energy 
cross section one has to consider an "off-mass-shell" amplitude. Since we are 
only interested in the cross section computed to an exponential accuracy, it 
does not matter what specific operator 0{(p) do we use to go off-shell. We 
choose 0{<f>) to be a low-power polynomial of (f). The analogue of eq. ( |l.lD in 
quantum mechanics would be 

a{E) = lmj dti,2e-^^(*^-*^) (O(0(ti))O(,^(t2))) • (2.2) 

One can decompose the Green function (...) as a sum over intermediate 
states with wave functions ^'„(0): 

2 



{om,))om2))) = E 



d(t>^o{4>)0{4>)Ki4>) 



■ 0{ti - t2)e-*^"(*i-*2) + e{t2 - ti)e-*^"(*2-ii)" , 
Integrating over ti^2 in eq. (|2.2| ) and taking the imaginary part we get 

2 



a{E) = TT5{En - E)V 



#vI/o(</')O(0)n(</') 



(2.3) 



(2.4) 



where V is the time volume of the process (actually to be cancelled by the 
omitted "flux factor" in eq. ( |2.2| )). We see thus that the quantum-mechanical 
"cross section" comes to a square of a matrix element between the ground state 
^0 and a highly excited one ^e- 

According to the Landau theory jl^ this matrix element can be calculated 
to the exponential accuracy as an exponent of the difference of two shortened 
actions, one with zero energy and the other with the given energy E: 



a{E) = exp2 



2E 



(2.5) 



where U{(j)) is the potential energy of the field. 

Eq. ( |2.5[ ) is directly applicable to the simple case of a one-minimum poten- 
tial; in this case one has to integrate in the first term from (j)min corresponding 
to the minimum of the potential, to 4> = oo, and in the second term from the 
energy-dependent turning point (pi where the integrand nullifies, to (j) = oo. 
For the square of the amplitude one has to double the result in the exponent. 
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We arrive thus to singular Euclidean trajectories: one starts at (Euclidean) 
time t = —oo from the minimum of the potential, goes according to the equa- 
tions of motion with zero energy to infinity (singularity) at some t = —T/2, 
then, starting from the same singularity, goes with fixed energy E to the cor- 
responding turning point (pi at t = 0. At the turning point = 0, and one 
can enter the Minkowski space to calculate the real-time evolution of the field. 
However, the Minkowski action is a pure phase, therefore to find the cross sec- 
tion one needs not know this part of the trajectory. Instead, one has to square 
the amplitude, i.e. to repeat the singular Euclidean trajectory in the opposite 
direction: starting from the above turning point (pi at t = 0, going with energy 
E to a singularity at t = +T/2, and ultimately returning to the minimum of 
the potential at t ^ +oo. The different branches of the trajectory are shown 
schematically in Fig. la. 

The Yang-Mills theory is however more like the double-well quantum me- 
chanics, therefore we have to generalize eq. ( p.5[ ) to the case of two-minimum 
potential as in eq. (2.1), allowing for the instanton transitions between the 



minima. To that end we have first of all to specify the initial and final states 
^o,£; entering eq. ( |2.4| ). In a double- well potential all levels are known to be 
split into symmetric (s) and antisymmetric (a) states. For levels deep inside 
the wells the corresponding wave functions can be written as superpositions of 
states localized in the left and in the right wells: 

^{^+1)±^{(P-1) 



V2 

^(0±l) = 5^iM4^. (2.6) 
V2 

To mimic in quantum mechanics the non-perturbative Chern-Simons chang- 
ing transitions of the Yang-Mills theory we choose the initial state with low 
energy (^'q in the notations of eq. (|2.4D ) to be localized near the left minimum 
and having approximately zero energy, 

vI/r"(<^) + ^~°(<A) 
^0(0) = ^ • (2.7) 

For the high-energy state we take a superposition which is predominantly 
localized in the right well and has energy approximately equal to E: 
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Of course at energies E higher than the top of the potential barrier one 
cannot say anymore that the state ( |2.8D is locahzed near the right minimum, 
but it is probably the best one can do to define continuously in energy the 
transitions between the left and right wells. In the Yang-Mills theory a much 
more clear trigger of the transition with the change of the Chern-Simons num- 
ber is given by the accompanying change of fermion chirality or baryon/lepton 
number violation. However, for energies less than the barrier top the state ^ e 
( |2.8D is indeed localized near the right minimum. For energies higher than the 
barrier matrix elements for the transition to any state are anyhow exponen- 
tially small, so the concrete definition of the final state does not make a great 
difference. 



With this definition of the initial and final states, eq. (2.5) has to be slightly 
modified. For energies higher than the potential barrier the modification is 
cosmetic: the singular trajectory described above has to start in the specific 
minimum (p = —1; it ends up in the same minimum - see Fig. la. For energies 
lower than the barrier the Minkowskian part of the trajectory starting from 
the turning point (pi (see Fig. lb) hits another turning point 03 where it again 
enters the forbidden zone and hence develops in imaginary (Euclidean) time. 
This branch (V) is a bounce resembling the kink plus anti-kink: at t = (chosen 
for symmetry reasons) the trajectory reaches the turning point 04 belonging 
rather to the right well. From this point one can proceed in Minkowski time 
observing the final-state field in the right well. However, the Minkowski action, 
being a pure phase, does not contribute to the cross section; in order to obtain 
it one has to repeat all the trajectory in the opposite order and to count the 
actions (with appropriate signs) along the branches I-V - see Fig. lb. We get 
for the actions along the different branches: 

J-oo 



TT TTT 1 r4>l=-\J^ + \/e I 

^ = ^ =2^1 ^ #v/(02-l)2-e, 

5^ = - / , W(02-l)2_e, (2.9) 

where we have introduced a dimensionless energy e = Sg^E. The final formula 
for the cross section is 

a{E) = exp [-S^ + S^^ - + S^^^ - S^^) . (2.10) 
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At energies higher than the barrier (e > 1) the branch V is absent, and there 
is no contribution of Sy to (7{E) so that we return to eq. ( p.SI ). 

Let us comment on a few features of eqs.(2.S. 2.1C| ). The actions to S^^ 
each diverge at (p = — oo however their sum is finite. To see that exphcitly we 
can rewrite their sum as 



I-IV 



+ 1) 



S 



= 5^ + 5 



IV 



S 



II 



S 



III 



+ 
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1)2 



> 0, 
(2.11) 

which apparently is convergent. At e — > these branches cancel altogether, and 
one is left with the piece which in this limit is an action on the infinitely 
separated instanton (kink) and anti-instanton (anti-kink), equal to A/Zg^. It 
gives the familiar Gamov suppression factor for tunneling at zero energy. We 
thus rewrite eq. ( 2.10| ) as 



(t{E) = exp 



(2.12) 



.S^-iv^E)- S^{E) . 

As the energy rises the behaviour of the cross section is determined by the 
interplay of two opposite trends. The tunnelling probability, represented by 
exp(— S^), increases with energy whereas the overlap of the initial low-energy 
state with the high-energy one, represented by exp(— S'^^''^^), decreases. This 



physical picture has been suggested several years ago in ref.[15|; now we are 
in a position to make it fully quantitative. At energies higher than the barrier 
there is no compensation for that trend from the side of , so the cross section 
will definitely decrease, and we thus expect a maximum of the cross section at 
energies of the order of the top of the barrier, e = 1. 



2.2 Calculation of the Cross Section 

The time dependence of the trajectories is given by 

t-to 



1 



(2.13) 



For e = (branches I and IV) we get trajectories going from —1 to — oo at 
t = ±T/2: 



b\t) = cthi(i + |), <^^^(t) = -cthi(t - ^) = <t>H-t)- (2.14) 
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At e ^ (branches II, III and V) the trajectories are given by ehiptic functions 
of the first and second kind. 

Actually in order to calculate the cross section ( |2.10| ) as a function of 
energy it is not necessary to know the trajectories explicitly: the shortenned 
actions contain less information. We first calculate the derivatives of shortenned 
actions in respect to e. At e < 1 we have (K(k) is the complete elliptic integral 
of the first kind 

^ de ^ de ^ de 



1 K(fc) K{k) / ^ 



At e > 1 we get: 



^ dlna _ _ ^dS^-^ _ 1 K(0 / 
^ de de ~ ~2 \' 

^ 0{e-'/% b(\^^= (2.16) 




8ei/4 V4' 2/ ' " V4' 2J r(f) 
Integrating these equations we get the cross section. At small energies we 
obtain a rising cross section, 



while at large energies we get a decreasing cross section: 

gHna{e) = -'—B (^-, - j +0(eV4). (2.I8) 

The maximum of the cross section is achieved at e = 1 corresponding to 
the top of the barrier, and the log of the cross section there is exactly twice 
less than at e — > where the suppression is solely due to the double instanton 
action: 

2 2 
g Inamax = --, e = I. (2.19) 
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2.3 "Square Root" Suppression 



The "square root" suppression of the cross sections at the maximum as com- 



pared to that at zero energies have been advocated by Zakharov |17] and Mag- 
It is 



giore and Shifman |18] from unitarity considerations (see also ref.|jT£ 
remarkable that we get it here in a different approach. Probably it means that 
the correct formulae respect unitarity. Moreover, we can show that this square 
root suppression is of a more general nature. Indeed, let us compare dS/de 
along the branches II, III and V for e < 1. According to a general theorem of 
the classical mechanics (see Appendix A) these derivatives are related to the 
periods of motion along these branches: 

dS'-'"" 



89' 



de 



u{4>) = (02 



-85- 



1)^ 



■^=4 
de 

--T1+T2. 



<^3 VUW 



(2.20) 



Here Ti,T2 are periods of motion along the branches V and 1 1 +11 1 respectively. 
We are going to prove that Ti = 2T2 for a wide class of potentials U{(f)). 

To that end let us consider the difference of integrals along two closed con- 
tours Fi and r2 in the complex (p plane (see Fig. 2): 



■I 







, _ , (2.21) 

VUicp) - e Jt, ^U{^) - e 

This difference is zero as far as the potential energy f7(0) has no singularities 
in the whole complex <j) plane, so that the integrand has only cuts at the turning 
points 01-4. According to eq. ( 2.20 ) the integral along the contour Fi is equal to 



Ti while that along F2 is equal to twice T2 , as there are two equal contributions 
from two cuts to that integral. Therefore we indeed get the relation 



dS^ dS'-'^ 



2, 



(2.22) 



de de 

Integrating this equation in e we arrive to the following important relation 



first noticed in ref. |2C]: 



S^{e) = 5^(0) - 2S^~^^{e). (2.23) 

Eqs.( 2!2^ , 2.23| ) are of a very general nature and are valid for any double- 
well potentials without singularities in the complex plane. (For the concrete 
potential considered eq. ( 2.22 ) follows directly from eq. ( p.l5| )). Moreover, we 
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think that eqs.(2.22,2.23) are valid also in field theories where tunneling may 
occur, Y-M and Y~M-H theories being an example. The analogue of eq. (2.21) 
in field theory would be the following contour integral in the complex plane of 
an arbitrary parameter s parametrising the solutions of eqs. of motion with a 
given field energy, (Pe{x,s): 



J (Px {d(j)E{x-, s) /ds) 



(2.24) 



^^d^x[\{V,<t^Ef + U[<l^]}-E 
One can choose the time t as the parameter s, but that is not necessary 



since (2.24) is re-parametrization-invariant. One can try to use this freedom 
to introduce a parametrization <j)E{x, s) such as to ensure that both the numer- 
ator and the denominator have no singularities in the complex s plane, except 
the cuts due to the turning points of the denominator. (In the one-degree-of- 
freedom case such parameter s is the "field" cj) itself). If that goal is achieved, 
one would prove the relations (2.22, 2.23| ) in field theory. We have not proven 
them in a general case, however we confirm these relations for the Y-M theory 
by a direct calculation in the next section - at least for small field energies. 
At e = 1 one reaches the top of the barrier where Sy vanishes. According 



to eq. (2.23) it means that 



S 



I-IV 



(1) 



5^(0) 



S 



inst 



1 2 

^3' 



(2.25) 



which reproduces eq. (|2.19| ) and gives the maximum of the cross section corre- 
sponding to the square root suppresion as compared to that at zero energy. 



Combining eqs. (2.23, 2.25) we can write: 



hia{e) 



_Sinst^ e = 1, 
-S'-'^{e), e>l. 



S'-'^'ie) > 0. 



(2.26) 



This equation shows a jump of the derivative at e = 1. It should be men- 
tioned though that in a parametrically narrow strip near e = 1 the semi-classical 
formulae are not valid anymore, and one would expect a smooth match between 



the two regimes of eq. ( 2.26 ). 
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2.4 Inst ant on Interactions 



Let us return to our concrete quantum mechanical example. To make contact 
with the previous work on the instanton-induced cross sections at low energies 
we recall that a conventional way to present the results is via the instanton- 
anti-instanton interaction potential, Umt ■ The cross section can be written as a 
Fourier transform of the time separation of the instanton and the anti-instanton 

[H, m, nil: 



a{E)= I dTexp [ET - Uint{T)] . 



(2.27) 



We define here Uint so that at large separations T it includes also twice the free 
instanton action. The integral over T is performed by saddle point method - 
that is how one gets the energy dependence of the cross section. Comparing 
eq. ( 2.27 ) with (t{E) found above (eq. ( 2.17] )) we obtain the instanton interac- 
tion potential at large separations: 



48Te-2'^ + 136e-2'r + 0{e'^'^] 



(2.28) 



The first term here is twice the instanton action, the second one repro- 
duces correctly the well-known kink-anti-kink leading-order attraction (see, 
e.g., ref.[^]). The third term is new: we can compare it only with the valley 
approach of Balitsky and Yung |2(:, 22, 27|. For the sum ansatz of a kink and 
anti-kink these authors get 



valley 
int 



9^ 



40e 



-2T 



+ 0{e 



-3T\ 



(2.29) 



which coincides with our result only in the leading order. To our mind, it means 
only that the sum of a kink and anti-kink is not too useful beyond the leading 
order. In an indirect way, however, our result for the next-to-leading kink- 
anti-kink interaction seems to be in accordance with another calculation, this 
time for the instanton-anti-instanton interaction in the Y-M theory. Indeed, as 
emphasized in refs.|26, 22, 27|, an instanton-anti-instanton confuguration can 
be obtained from a kink-anti-kink one through a chain of conformal and gauge 
transformations. The separation between the Y-M instantons R and their sizes 
pi^2 are related to the separation of a kink and anti-kink T as follows. 
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R' + pl + pl + JiR-' + pl + plY-^plpl 



(2.30) 



while the couphng should be replaced by the gauge coupling a according to 
the rule 



1 4 
?3 



47r 

a 



(2.31) 



Therefore, the kink-anti-kink interaction ( 2.28| ) can be transformed into 
the instanton-anti-instanton interaction. We get from eqs. ( ^ , |I30|) : 



Y-M 
int 



An 

{R) = - 
a 



qPipI I ^^PipUpi + P2) 



R^ 



(2.32) 



These terms are exactly those which are independently known today for the 
Y-M instanton interactions. It should be stressed that the last term has been 
computed from unitarity in a laborous two-loop calculation in the Y-M theory 
[ p^ ] 1^. It is remarkable that we reproduce the result of these hard calculations 
in such a simple way. 



3 Singular Yang— Mills Fields 

We are now turning to the Y-M theory and are going to describe the analogues 
of the branches 1-V (see Fig.l), assuming that the cross section of the 2 many 



processes are to the exponential accuracy given by eq. (2.10). A derivation of 



the generalization of the Landau formula to field theory using the first order 



formalism, can be found in the paper of lordansky and Pitaevsky |12]. We 
shall reproduce the well-known results for the instanton-induced cross sections 
at relatively low energies and get new results for high energies. In the case 
of QCD we shall get the cross section of isotropical multi-gluon production as 
function of energy: it decreases as approximately the ninth power of energy. 



^Recently the result has been confirmed in a non-less laborous calculation in ref.|28 ) for the case 
of the maximum-attraction orientation, though there is still a difference between refs. and 
for general orientations. 
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3.1 Exact Singular Solutions with Zero Energy 



In principle one could consider various types of field singularities but in this 
paper we restrict ourselves to the case which at low energies reproduces the 
instanton-induced processes. These processes are known to posses the 0(3) 
symmetry, meaning that the multiparticle production is spherically-symmetri- 
cal. One might well doubt whether it is natural for high-energy collisions 
and in a sense we prove that it is not, but at the moment we would like to 
follow the instanton-induced processes up to asymptotically high energies. For 
that reason we choose the singularity to be 0(3) symmetrical, meaning that at 
given times the field has a power- like singularity in r where r is the 

distance from the origin. 

For branches I and IV corresponding to E = Q exact solutions of the 
needed type have been already constructed by Khlebnikov |jll|: They are the 
usual BPST instanton (branch I) and antiinstanton (branch IV) in the singular 
Lorentz gauge with the scale parameter changed to —p^: 



T 
2 



t + ^ - p 



pi _ J.2 



{t + l-p 



+{t + 



{t + l-p) 



T 

t < , branch I, 

2 ' 



(3.1) 



2<, 



t-^ + p 



(t-l+p) r^+(t-^+p 



t > 



T 



branch IV, 



(3.2) 



where rj, rj are 't Hooft symbols 



Eq. ( |3.lD describes the evolution of the field with zero energy starting from 
zero field at t = — oo and getting to a singularity at i = — T/2. It is the 
analogue of eq. ( 2.14| ) in the quantum mechanical example. Eq. ( ^.2| ) describes 
the conjugate process: it goes from a singularity at t = r/2 to zero field at 
t = +00. In case of the EW theory with Higgses one can neglect the influence 
of the Higgses as far as the characteristic scale of the fields is p ^ where 



^We a grateful to C.Wetterich for a discussion of this point. 
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mw is the W boson mass. This is the case in two hmits we are now interested 
in: at E <^ rriw/a where we shall reproduce the well-known result, and at 
E ^ mw/a, see below. In the intermediate energy region one would have 



to solve the coupled Y-M-H system, otherwise eqs.(3J.,3^) are exact. At 



t = these solutions become singular at the origin, r = 0: 



4^Nr,-^)=^^4^, (3.3) 




4^^" f r, = ) " 'jf . (3.4) 



We notice that ( |3.3D and ( |3.4D coincide for spatial components Ai but 
differ in sign for the time component ^4. This is a gauge artifact, however. 
Indeed, one can perform a time-dependent gauge transformation 

A^, ^ [/t(r, t)^^?7(r, t) + iU'^d^JJ (3.5) 

eliminating the A^ components in both branches, I and IV. (In fact ^4=0 
is the adequate gauge to use: all results are physically more transparent in it; 
we started with a Lorentz gauge to make formulae more compact). It is easy 
to check that making a gauge transformation ( |3.5| ) we do not destroy the 
coincidence of the spatial components Ai at the singularity points. We cite the 
neccessary formulae in Appendix B. The time-dependent gauge transformations 
( |3.5D are however defined up to time- independent transformations. In order 
not to introduce new quantities, we use this gauge freedom to keep the Ai fields 
at t = ibT/2 exactly in the form given by eqs.(p^j3^): 



(3.6) 

In general, in the j44 = gauge the spatial components Aj(r,t) contain 
three 0(3) symmetrical structures, 



.0/ N \-A(r,t) .B(r,t) C(r,t) r _ 

A'^{r,t) = eijknj + {5a^ - nam)^^ + nani^^, n = -. (3.7) 
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In Appendix B we quote the relation of these structures to the field in the 
Lorentz gauge, and the Y-M action written in terms of the A, B, C functions. 

In order to find the fields along branches II and III corresponding to a 
non-zero field energy E field, one has to solve the Y-M equations (generally 
speaking, coupled to Higgses) in the time interval — T/2 < t < T/2 with the 
singular boundary conditions (|3.6| ) . As in the quantum mechanical example, 
the field energy E field is directly related to the time interval T. However difficult 
technically, solving the Y-M equations with given boundary conditions is a well- 
formulated problem which can be solved numerically. We solve it analytically in 
two limiting cases: (i) p <^T corresponding to low energies, aEp <C 1, and (ii) 
p ^ T corresponding to high energies, aEp ^ 1. In the first case we reproduce 
the usual instanton results. The second case corresponds to asymptotically 
large energies, and the results here are new. 



3.2 Reproducing Instanton Results 



In the case p <^T , i.e. low energies, we deal with a situation shown schemat- 
ically in Fig. lb. Therefore, one has first to find the singular solutions along 
branches II and III, then go to the Minkowski space and find the analogue of 
the turning point (p^. Then one has to construct again the Euclidean solution 
along branch V. Each time the boundary conditions for the solutions are deter- 
mined at the previous step, so that the procedure is straightforward. However, 
we have not solved the problem at arbitrary values of p/T but only in the lowest 
order in p/T. The explicit construction of the Y-M fields along the branches I 
- V is rather lengthy and we relegate it to Appendix C. We find the following 
expressions for the full actions along branches I+IV (with the actions along 
branches II and HI subtracted - we call it S^~^^ , see eq. ( p.ll| )) and along the 
branch V: 



qI-IV 



4tt 
a 



(3.8) 



'-'full — 



1 



In eqs.(3.8,OD we have introduced the time T2 = T — Ti which is the net time 
for branches II+III. 
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According to the general mechanical formula (see Appendix A), the energy 
of the field is the derivative of the full action in respect to the period of motion, 
which is T2 for branches II+III and Ti for branch V. Also, the field energy Efi^i^ 



should be the same on branches II, III and V. From eqs.(3^,3^) we have: 



E 



dS 



field 



I-IV 
full 



47r 

— 24-^ 

a r| 



"■'^full 



dT2 

— 24- -. 

a Ti - 5 



(3.10) 



dTi a [11 

We hence get a relation between the times along branches V and II+III: 



2T2 , where To 



-1/5 



In section 2 we have obtained the same relation for the case of quantum 
mechanics and presented arguments that it could be of a more general nature. 
Now we have demonstarted the validity of this relation "experimentally" in the 
Y-M theory, at least for large periods, or, equivalently, for small energies E. 

The correspondent shortenned action are (see Appendix A for the defini- 
tions): 



5, 



V 

short 



Sjull - EfieldTl 



An 



a 



V 967r J 



oI-IV 
short 



(3.12) 



(3.13) 



Adding these two actions in the exponent in accordance with eq. (2.12) we 
get finally the cross section as a function of energy at small values of aEp: 



a{E) = exp ■ 



47r 



a 



\ 967r J 



(3.14) 



This formula coincides with the well-known result for small energies [21 
pO| ] obtained for the instanton-induced total cross section in the saddle-point 
approximation. The one-loop correction to eq. ( p. 14 ) is also known |22, 31 
^ ]: it adds to the exponent of eq. ( |3.14| ) a term 



a V 967r 7 



(3.15) 



The 2-loop correction |2j] adds 



16 



47r 144 /aEpX 



8/5 



, 967r 

In +0(1) 



(3.16) 



We have not reproduced these terms directly by our new method (indirectly 
we got these terms from the quantum-mechanical calculation of the previous 
section by using the conformal-symmetry arguments). We would like to stress 
that what looks as a loop expansion in the conventional approach appears to 
be equivalent to finding singular solutions of the classical eqs. of motion. We 
believe that it is a well-formulated program which may be performed for all 
values of the dimensional energy aEp. 



3.3 Singular Solutions for High Energies 

Let us now consider the opposite case: T <^ p. It means that the time during 
which the field is developing is much less than the spatial spread of the field - 
see eq. ( |3.6D . During this short time the fields change only at distances r <^ p 
but cannot change significantly at r p. Therefore, the region r > p cancels 
out in the difference between the shortenned actions. For that reason one can 
neglect r as compared to p in eq. (|3.6D, so that the boundary conditions for the 



field simplifies. In terms of the A,B,C functions introduced in eq. (3.7) the 
boundary condition becomes: 

a(.,±|)=3. b(.,4)=c(,-.4)=-^. (3.17) 

We see that the structure B + C is singular and thus much larger than the 
structures A and B — C; therefore we shall neglect the second two and work 
with one structure B + C = D. The equations of motion for A,B — C and D 
will keep A and B — C negligible as compared to D as far as r will remain much 
less than p, which is the case of interest. Moreover, for the similar reasons one 
can neglect the spatial derivatives of -B + C and other less singular terms as 
compared to the terms B'^,C'^/r'^ in the action - see Appendix B, eq. ( [B.9D . 
The resulting action written in terms of the large structure D becomes: 



The corresponding equation of motion is 



D + —^ = Q (3.19) 
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with the boundary contitions 



D{r,±^) =-^. (3.20) 



The problem can be regarded as a simple one-degree-of-freedom mechanical 
one, with r viewed as a parameter of the trajectory. Eq. ( |3.19|) is integrable as 
there exists a conserved energy density: 

«(0 = ^ - Id'- (3.21) 

The time dependence of the field is determined by integrating eq. ( 3.21| ): 

fD{r,t) ^J) 

t= / (3.22) 

J-ip/r ^D^/Ar' - 2e(r) ^ ' 

Introducing dimensionless quantities, 



[8e(r)r2]i/4' r[8e(r)r2]i/4 ' 

we rewrite eq. ( |3.22 ) as 

f[8e(r)rY/^= /' ^^^^ (3-24) 
2r Jx \/x^ - 1 

Half-period of the motion, T/2, is found from eq. ( |3.24 ) as an integral from the 

initial configuration /3 up to the turning point being x = 1 in the new notations: 

f [8e(r)rY/^= f (3-25) 
It can be written as an equation relating the period and the energy density: 

dx _ pT 
Va;4 - 1 ~ r 

This equation has a solution for (3 and hence for e(r) for all values of the r.h.s 
The shortenned action density along the branches 11+111 is 



PI ^=,=-^- (3-26) 



64p^ 1 



,3 1 r(3 



^„ , dxVx'^ - 1. (3.27) 
p-^ Jl 

This action is apparently divergent, however one should not forget to sub- 
tract a similarly divergent piece corresponding to the zero-energy branches 1 
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and IV. To make the necessary subtraction we notice that zero energy corre- 
sponds to infinite period T and hence to /? — > cxd - see eq. (3.26). Putting 
/? = oo in eq. (3.27) we get the shortenned action density for branches I and 
IV: 



s'^''{r) = -^- (3.28) 

which is indeed independent of T, and should be subtracted from eq. ( 3.271 ). 
Of course, one could get eq. ( 3.2g| ) directly by finding the action density along 
the solutions ( p.l| , |3.2| ). To get the complete shortenned action standing in 
the exponent for the cross section we integrate the difference of ( 3.27] ) and 
( |3.28| ) over r. Using the relation ( p.26| ) we integrate over (3 instead of r, and 
obtain: 



4a Jo a \TJ vr Ji ^(3 \Ji 
4vr./i?(i i)\Vp\3/2 47r f p^l^ 



a 



We remind the reader that the above calculation is performed for the case 
p ^ T, hence eq. (|3.2S| ) describes a rapid fall-off of the cross section in this 
region. 

Integrating the energy density e(r) over r we get the field energy: 
P 3 47r /p\3/2 3 

dp ( ff^ dx ( fP dx P \ 



47r 3 (Biij)_y ^ 47r 05539 / p_y/' 



a T \T 



= ^^m ■ (3.30) 



One can immediatelly check that eqs.( p.29| , |3.3[l| ) satisfy the general relation 
between the shortenned action and the energy (see Appendix A): 

^ql-lV jp 

""^short _ field rp /o 01 \ 

dT dT ' ^ ' ' 
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The energy E field ( |3.30| ) has actuahy a double meaning. We have intro- 
duced it as the energy of a singular Euclidean field between the singularity 
points. It becomes purely potential energy at the turning point (t = 0) starting 
from where the field enters the classically allowed region and hence develops in 
real Minkowski time. Therefore, it is also the energy of the outgoing field or, 
equivalently, the energy of the final multiparticle state. 

Simultaneously, it is the energy of the few-particle initial state. To see 
that we notice that the cross section is defined via the imaginary part of the 
4-point Green function of the initial particles, see eq. (|]|). The high -energy 
asymptotics of a Green function in the background field is determined by the 
positions of the singularities of the field. In the case depicted in Fig. la there is 
one singularity at t = — T/2 for the amplitude and one singularity at t = +T/2 
for the conjugated amplitude. Therefore, the asymptotics of the 4-point Green 
function and hence of the cross section is 

exp [^i(-r/2) + E2{-T/2) - EiT/2) - E2T/2 - Sfuii{T)] 

= exp [-EpartideT - SfulliT)] (3.32) 

where Eparude = El + E2 is the energy of the initial 2-particle state and Sjuii 
is the full or Lagrange action of the field: 

Sfuu = sj^ii + Sf^ii — sy^ii — syjii- (3.33) 

On the other hand, the full actions are simply related to the shortenned 
actions provided the fields satisfy the equations of motion (see Appendix A) : 

SfuU = Sgfiori + Sl^^g^i — Sgl^j.^ — Sg[gj.^ — EfieidT = — EfieidT (3.34) 

(the quantity S^~^^ has been calculated in eq. (|]2|)). We get thus from 
eq. O): 



exp 



-EpartideT + + EfieldT 



(3.35) 



Integrating over the separation between the singularities T we obtain the 
energy conservation law, 

Epartide = E field = E (3.36) 



20 



whereas the cross section is given simply by 



cr{E) = exp 



-S 



I-IV 



(E) 



(3.37) 



where the shortenned action S^~^^ should be expressed through the energy. It 
should be mentioned that in a case depicted in Fig. lb there are additional field 
singularities at t = ibTi/2 and a similar analysis is slightly more complicated. 

Expressing the time between the singularities T through the field energy 
( |3.3[1| ) and substituting it into eq. ( 3.2g| ) we get finally 



cr{E) = exp 



1.316 1 ^ ' 



(3.38) 



a \ An J 

which is valid for aEp S> 1 and describes a decreasing cross section at large 
values of this parameter. 

The result of this subsection can be formulated also in terms of the instan- 
ton-anti-instanton interaction potential, Uint ■ We find at small separations T 
a strong repulsion. 



s 



I-IV 



ET 



—0.3693 ^ 

a \T 



3/2 



T ^p, (3.39) 



as contrasted to the usual attraction at large separations, see eq. ( 2.32 ). We 
believe that it is a reasonable result, if one defines Uint not as the action of 
an arbitrarily chosen Euclidean configuration like the valley (we have shown 
in section 2 that the valley leads to the wrong Umt already in the next-to- 
leading order), but through unitarity, as in eq. ( ^^ ). Physically, the effect ive 
repulsion of instantons at small separations appears here as a result of a rapidly 
decreasing overlap between low-energy and high-energy states. From the purely 
Euclidean viewpoint a closely situated instanton and anti-intstanton resemble 
a vacuum state and therefore have to be effectively strongly repulsive if we wish 
to read off a non-perturbative contribution from that configuration. Eq. (3.39) 
is good news for the QCD instanton-vacuum builders. 



4 Energy Behaviour of the Cross Section 

We have established the behaviour of the cross section for multi-gauge-boson 
production as a function of a dimensionless quantity aEp where p is the spatial 
size of the Y-M field, at small (eq. ( p.l4D ) and large (eq. ( |3.38| )) values of 
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this parameter. Applications of this result are different in the QCD case where 
scale invariance is preserved at the classical level, and in the electroweak theory 
where it is broken from the beginning. We discuss the two theories in turn. 



4.1 QCD 

In QCD the scale p is arbitrary, and one has to integrate over it. Naturally, 
p will be then fixed by the maximum of the cross section as function of the 
dimensionless parameter aEp. We know the behaviour of the \oga{E) at small 
and large values of this parameter, and we plot it in Fig. 3, where we have added 
the next-to-leading term ( 3.15| ) at small aE'p. The low-energy and high-energy 



curves intersect and give the maximum at a point which is remarkably close to 
a point marked by a circle in Fig. 3, which is the position of the maximum that 
we would expect independently of the above calculations. 

Indeed, in Section 2 we have presented general arguments in favour of the 
relation ( |2.23| ) , and in Section 3 we have confirmed it at low energies for 
the Y-M case. If we assume eq. ( 2.23| ) to be correct at all energies up to the 



sphaleron top of the barrier, we would get that at the maximum corresponding 
to the top of the barrier the cross section is exactly the square root of the cross 
section at zero energy. It means that the low- and high-energy branches of the 
log of the cross section should intersect at 

\ogamax = -• (4-1) 

a I 

To what value of the dimensionless parameter aEp does this maximum cor- 
respond? Our arguments here are even more shaky, however, one could specu- 
late that the maximum occurs at energy equal to the mass of the sphaleron. In 
the massless Y-M theory we are now dealing with, the sphaleron mass can be 
estimated as the potential energy of the instanton field exactly in the middle 
of the transition when the instanton passes the Nqs = 1/2 point. We have: 



M,^h[p) = A^vr / drr^ 



4g2 J ^ (^2 + ^2 + ^2)4 



(4.2) 

a pl6 ^ ^ 



t=o 

We find then that the maximum of the cross section should be achieved at 
the value of the dimensionless parameter 

(^) = A = 0.1875. (4.3) 
V 47r J^ax 16 
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This is the point marked by a circle in Fig. 3, and we see that it is remarkably 
close to the intersection of the low- and high-energy curves. Higher order 
corrections to both curves should probably move the intersection point exactly 
to the position calculated here. 

The fact that the spatial size of the classical field p scales down as 1/E at 
the maximum means that quantum corrections to the classical calculations of 
this paper should be controlable and small at — > oo. The one- loop quantum 
corrections, as usually, should make the gauge coupling run. Also, one would 
expect a large prefactor owing to the zero modes about the classical trajectory 
- something like (27r/a)'*^'', where Nc is the number of colours. The prefactor 
a starts to "run" only at the 2-loop level. Combining renormalization- group ar- 
guments with the result ( [4.1[ ) we expect the isotropical multi-gluon production 
cross section to be 

-...(^) = 1^ (I)' ("4)"'" [' + ° (mITX))] 

where h = llA''c/3 — 2Nf/2> ~ 7 is the 1-loop Gell-Mann-Low coefficient and 
A is Aqcd ~ 200MeV. The coefficient c and the power p are determined 
by the second loop Gell-Mann-Low coefficient and by the Green functions of 
the initial particles of the process; that is where the non-universality of cross 
sections comes in. A calculation of these constants seems to be feasible. 

We thus arrive to a prediction that the total cross section for isotropi- 
cal multi-gluon production in QCD should fall off as approximately the ninth 
power of their aggregate energy. However a presumably large prefactor and the 
peculiarity of the events might help to make such processes observable; this 
subject deserves a more thorough study. 

It is worth mentioning that the momentum distribution of the produced 
multi-gluon state can be easily found by solving the Y-M eqs. in Minkowski 
time starting from the "sphaleron" field ( [4. 21 ) . Since that field has only 
one parameter p related to energy through eq. ( ^.S]) , the gluon momentum 
distrubution will have a scaling behaviour with energy, modulo logarithmic 
corrections. 

4.2 EW Theory 

In this case the size of the classical field p is determined not only by energy 
but also by the Higgs v.e.v. v. At low energies one has to multiply the aEp- 
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dependent cross section ( [3.14D (with the addition of 1-loop ( 3.15| ) and 2-loop 
( 3.16| ) corrections) by the factor 



exp(-27rV2^2)_ (45) 

Integrating over p one gets the familiar expansion of the total baryon number 
violating (BNV) cross section in terms of the dimensionless energy measured 
in units of the sphaleron mass |24|: 



(^tot{E) = exp 



4vr , 
—F{e) 
a 



E 



F{e) = -1 + - -e^l^ - -6^/3 [ine + 0(1)] . (4.6) 

8 16 16 

The growth of the function -F(e) has lead in the past to hopes that the BNV 

processes may become unsuppressed at energies of the order of the sphaleron 



mass, m\Y /oi, corresponding to e ~ 1. However, unitarity arguments of refs.|17, 
p^ , indicated that the BNV cross section induced by instantons should not 
rise beyond the square root of the cross section at zero energies. We now arrive 
to the same conclusion from our singular solutions. 

At very high energies {E ^> ra\Y /a) the size of the field p is determined by 
the factor ( [3.38 ) rather than by the Higgs factor like (^]^) , meaning that 
p ~ \/aE <^ l/m\\r. Therefore one can neglect the Higgs influence and arrive 
to the conclusions of the previous subsection, viz. that the cross section does 
not rise above the value of exp(— 27r/a) ~ 10~®^. 

At energies of the order of the sphaleron mass {E ^ rriw/a) one cannot 
neglect the Higgs effects and has to find singular sloutions of the Y-M-H equa- 
tions along the branches I-V. If our general arguments based on analyticity 
are correct, we would again expect that the maximum occurs at exactly the 
sphaleron energy and corresponds to the square root suppression. 

The momentum distribution of the produced W and H bosons should follow 
from the Minkowski development of the sphaleron - the "fall" of the sphaleron. 



This distribution has been studied numerically in refs. |34, 35 1. We argue now 
that the same distribution would happen at high-energy collisions, however its 
probability would be of the order of 10~®^. 
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5 Massless Xcp^ Theory 



For future reference we would like to cite here the results for the massless Xcj)"^ 
theory whose Euclidean action is 



S = Jd'x(^^{d,cPf + ^^^^ (5.1) 



Instead of the singular instanton one has to start here with the singular 
Lipaton which gives an 0(3) symmetrical zero-energy field singular at t = 

±T/2 g-. 

Since there is no tunneling and consequently no branch V in this theory, one 
gets a monotonically decreasing cross section as function of the dimensionless 
quantity XEp. Probably it is in accordance with the theory being not asymp- 
totically free. Calculations of Section 3 can be repeated without any serious 
change (in fact they are much more simple). For small XEp the cross section 
has been actually calculated in ref.|jll|], and we confirm it: 



For large XEp the leading term can be obtained immediately from eq. (3.38) 
by changing notations. Indeed, introducing a new variable D(r,t), 

<Pir,t) = ^^, (5.4) 



2A 



we get for the action (|5.lD : 



which coincides with ( p.l8|) up to the overall factor. The boundary conditions 
for D(r,t) as imposed by eq. (5.2) also coincide with those of the Y-M case 



( |3.2[1| ) . Hence the behaviour of the cross section can be immediatelly obtained 
from eq. ( |3.38| ) by a suitable change of constants. We get: 

167r2 f XEp\^/^ , , , 

\ogatot = ^1.119 f^j , AEp»l. (5.6) 
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Eqs.( ^.3|j5.6| ) exhibit a monotonically decreasing cross section as function 
of XEp; integrating over p one falls into the trivial maximum at XEp = 
corresponding to zero particle production, which is naturally more probable 
than multi-particle production. It is a prerogative of field theories with non- 
trivial topology to have a non-trivial behaviour of classical cross sections even 
in the massless limit. 



6 Conclusions 

We have elaborated the use of singular Euclidean solutions to describe semi- 
classical multi-particle production in field theories. We have shown that using 
instanton-like 0(3) symmetrical singular solutions one reproduces the usual 
formulae for the instanton- induced cross sections for low energies. However the 
present approach allows to solve the problem at all energies. We have found 
analytically a new family of singular solutions labelled by p/T, which govern 
the high-energy behaviour of multi-particle production - both for the Y-M and 
A(/>^ theories. We have presented new arguments (not based on unitarity but 
rather on analyticity) that the maximum of the instanton-induced cross sections 
in QCD and EW theories is still exponentially small in the coupling constant, 
with the coefficient corresponding exactly to the "square root" suppression 
as compared to zero energy. An interpolation of the low- and high-energy 
branches of the cross section give numerically this "square root" point to a 
surprisingly good accuracy - see Fig.3. Simultaneously our result means that 
there is effectively a strong repulsion between instantons at small separations. 
As a by-product of our study we get a non-perturbative formula for the 



isotropic multi-gluon production in QCD, eq. (4.4). Presently unknown coef- 
ficients and powers of the logarithms in that formula seem to be calculable, 
and deserve further study. It would be also extremely useful, if possible, to 
understand eq. (|4.4|) in terms of some clever summation of Feynman graphs for 



a gluon branching process - we mean the approach initiated by Cornwall |36] 
and Goldberg Q. 

As to the applications to the baryon number violating processes in the EW 
theory, our result is in the line with other people who in the recent years ar- 
gued from various angles that the BNV instanton-induced cross section has to 
be small at all energies. However we feel that it is not the end of the story 
but, rather, its beginning. It should not be surprising that the isotropic multi- 
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particle production is suppressed at high energies. One should expect that at 
high energies the events have at most the axial symmetry. Therefore, one should 
look for solutions with singularities of a different kind - not just 0(3) symmet- 
rical. The same applies to other cases where we expect that multi-particle 
production might be described semi-classically, - such as multi-jet production 
and low-x structure functions in QCD, multi-pion production in heavy-ion col- 
lisions, etc. 

Our final remark concerns an alternative approach to semi-classical multi- 
particle production [^8| in which the few initial high-energy particles are used 
as a (weak) source for classical fields. How can this approach be reconciled 
with our? We think that including the source explicitly, one has to look for 
anomalous solutions of field equations, which do not possess singularities what- 
ever small is the source but which develop a singularity when the source is 
put identically to zero. We have observed such scenario to happen in quantum 
mechanics and a d=2 massive a model, and believe that it is a general case. 
However, in a d=4 Yang-Mills theory the hopefully equivalent approach based 
on singular Euclidean solutions is seemingly more simple. 
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Appendix A. General Relations for the Ac- 
tion 

Let S full be the Lagrange action: 




(A.l) 



Using the energy conservation, 




(A.2) 
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which is fulfihed if the field (p{t) satisfies the eqs. of motion, we can present 
Sfuii as 

Sfuli = ^J dt^^ + Ejdt = 
^ J dcP^U{<P) - e + ET = Sshort + ET (A.3) 



where T is the period of the motion. Eq. ( A.3 ) introduces the so-called short- 
enned action, 

Sshort = j pdq. (A.4) 

In the main text we make intensive use of the general identities following 
from eq. (|A.3| ): 

dSshort m dS short rridE / a r\ 

and 

^ = E. (A.6) 

Branches I and IV correspond to zero energy, so that the full actions along 
these branches coincide with the shortenned ones: 



Sfull — Sfull — Sghort — ^Ihorf {■^■'^) 



The relations ([A.5|, A.6) apply directly to branch V in which case the period 



is denoted by Ti, see Fig. lb: 

dS'^hort _ dSYhort _rj.dE dSj^ii 

The period of motion along branches II+III is denoted by T2 (it becomes 
equal to the full time between the singularities T for energies above the barrier), 
and we have: 

S'/ul/'' = SlL'/' + ET, (A.9) 

We usually combine this action with that along branches I and IV in order 
to cancel the divergencies. Denoting 
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qI-iv — ni+iv qII+iii qI-iv — qI+iv qII+iii {^^n\ 

^full — '^full '^full ' — short short ' l^A.iUj 

we have 

Sj-il"" = S'-'"" - ET2 (A.ll) 

so that 

dS'-'^ dS'-'^ dE dS]-;^,'' 

Appendix B. Gauge transformation to the 
A4 = Gauge 

Let us consider an instanton-like field in the Lorentz gauge: 

^m(^) = 2<.a:.^(a^') (B.l) 



or 



A^ = At^ = {nT)r^r^ + f), 
Ai = Af-^ = [-Tit + i ((nr)ri - m) r] ^(r^ + t^). (b.2) 



We make a hedgehog gauge transformation 



A!^ = C/t(r, t)A^U{r, t) + iU'^d^U, U = exp [i(nr)P(r, t)] (B.3) 

in order to ehminate the A'^ component. The profile function of the gauge 
transformation is found from the equation 

BP 

Ai = iUdiU^ = (nr)— = (nr)r$(r2 + 1^) (B.4) 

which can be integrated once the function <I>(r^ + 1^) is given. 

With P(r, t) known, the gauge-transformed spatial components A'^ take the 
0(3) symmetrical form, 



l-A{r,t) ^B{r,t) C(r,t) r ,^ ^, 

4«(r,i) = eijkUj ^^^ + {5ai-nani)^-^+nani^-^, n = - (B.5) 
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where 



A{r, t) = cos 2P + (cos 2Pr + sin 2Pt)2r$, (B.6) 



B{r, t) = - sin 2P + (sin 2Pr - cos 2Pt)2r^, (B.7) 



C(r,t) = -2r 



or 



(Bi 



If the 't Hooft symbol r] instead of f] is used in the definition of the field 
( p.lD , one has to change t — > — t in all the above equations. 
The Y-M action rewritten in terms of these structures is 



Sy.M = ^,jd'x{F^^:f = \^jdtjdr 



2 



^ ./2 , 0/2 , {A^ + B^-lf , 2C{A'B-AB') , C^A^ + B^) 
+ ^ +B + + + 



(B.9) 



Appendix C. Singular Y— M fields at Low 
Energies 

C.l Branches II and III 

At low energies we do not expect large deviation on branches II and III from the 
zero-energy solutions ( |3.1| , |3.2D which are just singular instantons. To save the 
space we introduce the following notations. Let j4^"**(t + T/2 — p) be a singular 
instanton in the Lorentz gauge defined by eq. (3.1). Similarly, ^J^"**(t — T/2 + /9) 



is a short-hand notation for the singular anti-intstanton given by eq. ( |3.2| ). The 
time arguments of the fields are what we wish to stress here. We next gauge 
transform the fields to the ^4 = gauge according to the formulae of the 
Appendix B. We denote the resulting spatial components as Ai{t + T /2 — p) 
and Ai{—t + T/2 — p) (note that we have changed the sign of the argument 
for the case of the anti-instanton field in agreement with the procedure of the 
Appendix B). 

We next notice that the net time for branches II and III is T2 = T — Ti 
(see Fig. lb). The field along branch II has to stop at a turning point at time 
—Ti/2, and the field along branch III has to stop at time Ti/2. Since these are 
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the same turning points, we can temporaly identify them with a single point 
t = 0; the singularities occur then at time zizT2/2. 

To get the fields satisfying the boundary conditions at the singularity points 
we introduce the field 



A,it + T2/2 + p), t<0 
" ^ Ai{-t + T2/2 + p), t>0. 

The two branches match at t = 0, however their time derivatives there have 



opposite signs. Eq. ( C.l ) satisfies the Y-M equation of motion everywhere 
except the point t = 0: 

= 2MT2/2 + p)5{t) = 5{t)0{p'M). (C.2) 

Therefore, we have to modify in the vicinity of the turning point at t = 0. 
We look for the solutions along the branches II and III in the form 

A{t) = A'i{t)+B,{t) (C.3) 

where the additional field Bi satisfies the zero boundary conditions at the sin- 
gularity points, ±T2/2, and the equation of motion: 

° = ^-^^^^ = ^ + i-d! + ICiA^))B + OiB^) (C.4) 

where — 9| + /C is the quadratic form of the Y-M action. 

To the accuracy we are now interested in, the field Bi is of the order of 



p^/T^. Therefore, we can neglect the non-linear terms in eq. ( C.4 ) and the 
shifts by p in the arguments of the fields. Furthermore, in the region of i « 
one can put A^ = in the quadratic form fC, since A^ is also of the order 
of pV^I in this region. Solving eq. ( |C.4D at t ~ and taking into account 



eq. (|C.2| ) we get: 



Bi{t^0) = Ai{\t\+T2/2). (C.5) 

The resulting field A^ -|- Bi is smooth at t = 0; it is shown schematically in 
Fig. 4a. Eq. ( |C.5| ) gives the only solution of the equation of motion to the needed 
accuracy, which is even in t and decreases in the direction of the singularity 
points. One could, in principle, add to eq. ( p.5| ) a solution of the homogeneous 
equation. Such solutions would be the zero modes of the singular instanton; 
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their addition would correspond to a shift or a distortion of the singularity 
points and therefore would not satisfy the necessary boundary conditions at 
±T2/2. 

We thus obtain the field at the turning point t = 0: 

Al''''\r,0) = 0, Al''''\r,0) = 2Mr,T,/2), (C.6) 

Let us now calculate the full action along the branches II and III, subtracting 
those along the zero-energy branches I and IV: 



'^'''^ (C.7) 

Using equations of motion and gauge invariance this action can be calculated 
without undue difficulty. Let us add and subtract the action computed on 



(C.l) in the time interval between =pT2/2. We write: 



Wi = 2 f ^^^^ dt f d^xC[Ai{t + T2/2)] 

J —CO J 

-2 r dt f d^x£[Ai{t + T2/2)], 
J-T2/2 J 

W2 = -J ^ dt J d^x^ClA'^ + B]-C[A^]y (C.8) 

In Wi we change the time variables and use the fact that the Lagrange 
density is even in time. We obtain: 

Wi = 2 dt d^xC[A,{t)]. (C.9) 
JT2/2 J 

This integral can be calculated in two ways: (i) using the equation of motion 
for the Ai field and (ii) using the gauge invariance. The first method leads to 
the equation: 

Wi = -^l d^xA^{r,T2/2)A'^{r,T2/2). (C.IO) 
Using the gauge invariance we can return to the Lorentz gauge and obtain: 
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In W2 we expand the integrand in Bi neglecting terms 0{B^) and 0{B'^) 
and use the equations of motion for and Bi (eqs.( pr2 , C.4| )). We have: 

1 rT2/2 



W2 = -^ r dt I d^xBf{r,t)A^{r,T2)6{t) 

= j d^xA'}{r,T2/2)A'^{T,T2/2) = Wi (C.12) 
where eq. ( ClOj ) has been used in the last line. Adding up the two pieces 



we thus obtain the full action along branches II+III (with branches I and IV 
subtracted) : 



(C.13) 



C.2 Branch V 



We now go back to the turning point at t = where we have already found 
the field - it is given by eq. ( |C.6D . What we called t = corresponds actually 
to two points: t = ibTi/2. At these points the field enters the allowed region 
and develops in Minkowski time. The solution of the Minkowski equations of 
motion is 

^Afinfc ^ Ai{T2/2 + it) + Ai{T2/2 - it). (C.14) 

Indeed, this is a solution of Minkowski eqs. of motion for not too large t since Ai 
itself is a solution of Euclidean ones, and ( C.14| ) is evidently obeying boundary 
conditions (^ 



The field ( C.14 ) describes the outgoing gluon field which corresponds to 
multi-gluon production without changing of the Chern-Simons number. In 
order to change it we need one more tunneling transition corresponding to 
branch V of Fig. lb. Normally the Minkowski solution develops between the 
turning points 0i and ^3 (as it was in our quantum mechanical example), and 
determines the initial field configuration (f)^ for branch V. However this is not the 
case in the lowest approximation we are dealing with. In this approximation 
the two turning points, ^3 and 0i, coincide, so that branch V starts at the 
same field where branch II ends. It should be also mentioned that far from the 
centers in the leading approximation there is no difference between the field Ai 
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of the singular and of the non-singular instantons. We shall therefore denote 
the instanton field in the = gauge by Ai as before. 

We have thus to find branch V as a solution of the Euclidean eqs. of motion 
with the boundary conditions given by eq. ( |C.6| ): 

AY{±Ti/2) = 0, AY{±Ti/2) = 2Ai{T2/2). (C.15) 
Again we are looking for a solution in the form 

AY = A'^i+Bi (C.16) 

where this time A^ is the field of the usual instanton in the ^4 = gauge at 
t < and that of an anti-instanton at t > 0; the centers of both are shifted by 
AT, to be found below: 

^0 = / Mt + Ti/2-AT), t<0 
" \ Ai{-t + T2/2-AT), t>0. 

(see Fig.4b). 

The shift in the positions of the instantons, AT, should be determined 
from the region near — Ti/2. In this region the field Bi obeys linearized eqs. of 
motion following from the quadratic form of the action in which one can neglect 
the instanton field being small in this region. The boundary conditions for Bi 
follow directly from eqs.( p.l5 . C.17| ): 



Biit = -Ti/2) = 2AiiT2/2) - Ai{AT), Bi{t = -Ti/2) = -ii(AT). (C.18) 

The field Bi should not increase in the direction to the center of instantons 
- then it does not include zero modes corresponding to the distortion of the 
instantons. We can fulfill this requirement only if AT = T2/2. Indeed, in 
this case the solution of eq. of motion for Bi with the boundary conditions 

(ED 



IS 



Bi = Ai{-t + {Ti-T2)/2), (C.19) 

which clearly decreases for increasing t and hence does not contain the zero 
modes which would lead to the shift of the instanton position. 
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Let us turn now to the region near t = 0. We cannot apply here immediately 
the method used above for branches II and III since the field A'^ is not small 
now. However it is gauge equivalent to a small field A^: 



A^iit) = U^r)A^U{r)+iU^diU 



(C.20) 



where U{r) is a gauge transformation with the winding number equal to unity. 
One can apply the same considerations as above to the gauge transformed fields 
A^ and Bi = UBiU^. As a result we obtain the following expression for Bi near 
t = (cf. eq. (|Cl|)): 



Bi = U\r)Ai (\t\ + 



(C.21) 



The solution of Minkowski eqs. of motion, which starts at the turning point 
at t = 0, is (cf. eq. dCl^ )): 



A 



Mink 



(t) = C/t(r) ^Ai i^it + + [-it + ^^jj Uir) + iU^diU. 

(C.22) 

It describes gluon production around the minimum with Chern-Simons number 
equal to one. 

Now we can calculate the action along the branch V. The full action, 



'-'full — ~2 , 

9 J-Ti/2 



f [ d\C[A^, + Bi], 

J-Ti/2 J 

can be expressed as a sum of twice the instanton action. 



S^^ = — I I SxC[Ai 



(C.23) 



(C.24) 



g J-oo J a 
and corrections of four types. The first is due to the fact that we have to calcu- 
late the action not along the whole t axis but in the time interval (— oo, — Ti/2): 



2 /-^i/s , , 
A5i = — 2 / ' 

g J -co 



Ai[t + 



T1-T2 



a \T2, 



(C.25) 



The second is a similar defect of the action but coming from the time interval 
(0, 00): 



AS: 



2 ~ — 2 , 
5 -'0 



Ai t + 



T1-T2 



An 



a 



T1-T2 



(C.26) 
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The third and fourth corrections are the contributions due to the Bi field: 

ASs + AS4 = \ I d'x + Bi] - CiA'i] I . (C.27) 

9 i-Ti/2 J ^ ' 

As above, we shall use here the eqs. of motion both for the field and 
Bi, and get two contributions. One is the contribution of the boundaries at 
t = ±Ti/2: 



= ^ (Fx 



B,{t)A, ( t + ^^^) + lBi{t)B,{t) 



The other is the contribution of the 6 function at t = 0: 



t=-Ti/2 



ASi. 



(C.28) 



A54 = l/ 



(fx 



Bi{t)Ai ( t + ^^^) + \Biit)Bi{t) 



AS2. (C.29) 



t=o 



Adding up the four corrections we obtain finally the full action along branch 



V: 



S 



full 



47r 



a 



1-6 



T2 



T1-T2 



(C.30) 
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Figure Captions 



Fig.l. Singular trajectories below (a) and above (b) the sphaleron 

Fig. 2. Contour integration proving the relation between the periods, eq. ( |2.22| ) 

Fig. 3. The behaviour of the log of the cross section as function of aEp at 
small (dashed curve) and large (dotted curve) values of this parameter. The 
circle shows the expected maximum of the cross section 

Fig. 4. A schematic view of the singular solutions along branches II and III 
(a) and the bounce solution along branch V (b) 
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